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Summary. A semi-classical mechanism which leads to the isotropization of the Mixmaster 
Universe is developed. A wave function of this Universe, which has a meaningful probabilistic 
interpretation, is constructed and it describes the evolution of the anisotropics of the Universe 
with respect to the isotropic scale factor, which plays the external observer-like role. We show 
that, once large volume regions are investigated, the closed Friedmann-Robertson- Walker 
configuration is deeply privileged. 

Quantum cosmology denotes the application of the quantum theory to the entire Universe 
[1]. It can be then viewed as a natural arena to investigate as part of a more general drive 
to understood quantum gravity. In canonical quantum gravity, the quantum state of the 
system is generally represented by a wave functional describing the dynamics of the three 
metric hij as well as matter fields 4>, i-e. •f' = ^[hij{x), 4>{x)\. This state is defined on the the 
Wheeler superspace and, since of the diffeomorphisms invariance of GR, it has no explicit 
dependence on time. In particular, it (formally) satisfies the Dirac quantum implementation 
of the first-class constraints of GR. 

A complete quantum theory of gravity is not yet available and therefore this problem is 
not properly defined. To overcame such a feature the fields are usually restricted (by hand) 
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to a finite dimensional subspace of the superspace, i.e. we deal with the minisuperspace 
representation. Quantum cosmology is explicitly defined as the minisuperspace quantization 
of homogeneous (finite degrees of freedom) cosmological models. From the peculiarity of 
the system-Universe, fundamental interpreting difficulties of the wave function of the Uni- 
verse •f' arise. The question about the interpretation (i.e. extracting physical statements) 
of quantum cosmology clearly appears as soon as the differences with respect to ordinary 
quantum mechanics are addressed [2, 3]. The standard interpretation of quantum mechanics 
(the Copenhagen one) involves the following basic assumptions, (i) There exist an external 
observer to the quantum system, i.e. the model under investigation is not genuinely closed, 
(ii) Predictions are probabilistic in nature and performed by a measurement of an exter- 
nal agency, (iii) Time plays a central and peculiar role. By contrast quantum cosmology is 
defined by the following features, (i) The analyzed model is the Universe as whole, i.e. it 
is closed without external observers, (ii) No external measurement crutch is available and 
an internal one can not plays the observer-like role since of the extreme conditions a very 
early Universe is subjected on. (iii) The time coordinate is not an observable in GR and at 
quantum level it is known as the problem of time. 

The most developed idea to solve these features relies in accepting that a meaningfully 
interpretation of the wave function of the Universe can be only formulated at semi-classical 
level. More precisely, it is only possible to quantum-mechanically interpret a small subsystem 
of the entire Universe, i.e. in the domain where at least some of the minisuperspace variables 
are semi-classical in the sense of the Wentzel-Kramers-Brillouin (WKB) approximation. 

In this work such a scheme is implemented to the most general homogeneous cosmologi- 
cal model, i.e. the Mixmaster Universe. The relevance of this model relies on the fact that a 
generic solution of the Einstein equations toward the cosmological singularity is formulated 
by a collection of causal independent Bianchi IX horizons [4] . In particular, a wave function 
of the Universe which has a clear probabilistic interpretation when the isotropic scale factor 
a of the Universe is regarded as semi-classical is obtained. It describes the quantum evolution 
of the Mixmaster anisotropies and its dynamics is traced with respect to o, which can be re- 
garded as a semi-classical variable as soon as the Universe expands enough. The main result 
is that the wave function of the Universe is spread over all values of anisotropy near the cos- 
mological singularity but, when the radius of the Universe grows, it is asymptotically peaked 
around the isotropic configuration. The closed Priedmann-Robertson- Walker (FRW) cosmo- 
logical model is thus the naturally privileged state fax enough from the classical singularity. 
A semi-classical isotropization mechanism for the Mixmaster Universe is then predicted. 

The dynamics of the homogeneous cosmological models (the Bianchi Universes) is sum- 
marized, in the canonical formalism, by the scalar constraint (for reviews see [5]) 



where the potential term V{l3±) accounts for the spatial curvature of the model. Here k = 
8nG, the variable a — a{i) describes the isotropic expansion of the Universe and its shape 
changes (anisotropies) are associated to f5± = p±{t). The phase space of this model is thus six 
dimensional and the cosmological singularity appears for a — > 0. In the Universe dynamics we 
have assumed the matter terms to be negligible with respect to the cosmological constant A, 
i.e. the isotropic potential U (a) reads U(a) = —a/AK + Aa'' / k. As matter of fact, far enough 
from the singularity, the cosmological constant term dominates on the other ordinary matter 
fields and such a contribution is necessary in order to the inflationary scenario takes place 



As we said, a correct definition of probability in quantum cosmology can be formulated by 
distinguishing between semi-classical and quantum variables [2] . More precisely, the variables 
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[6, 7]. 
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which satisfy the Hamilton-Jacobi equation are regarded as semi-cleissical and is assumed that 

the quantum variables do not affect the dynamics generated by the semi-classical ones. In 
this respect we claim that the quantum variables describe a small subsystem of the Universe 
and is then natural to regard the isotropic expansion variable a as the semi-classical one 
while considering the anisotropy coordinates /3± (the two physical degrees of freedom of the 
Universe) as the purely quantum variables. We are thus requiring ab initio that the radius 
of the Universe is of different nature with respect to the anisotropics. To implement such a 
picture, the wave function of the Universe If' — 'I^{a, j3±) is assumed to be [2] 

1^ = 1^oX = A{a)e^''^'^\{a,l3±). (2) 

This wave function is WKB-like in the a coordinate and the additional function x depends 
on the quantum variables f5± and only parametrically, in the sense of the Born-Oppenheimer 

approximation, on a. 

The canonical quantization of this model is achieved by the use of the Dirac prescrip- 
tion for quantizing constrained systems [8J, i.e. imposing that the physical states are those 
annihilated by the self-adjoint operator 7^ corresponding to the classical counterpart (1). 
Considering (2), we obtain from the quantum operator version of (1) the Hamilton-Jacobi 
equation for S and the continuity equation for the amplitude A 

-KA{S'f + aUA + Vc = 0, i [a'S') ' = 0, (3) 

respectively. Here the prime denotes differentiation with respect to the scale factor a and 

Vq = kA" is the so-called quantum potential, which in this model is negligible far from 
the classical singularity even if the ft ^ limit is not taken into account (see below). As 
usual S{a) defines a congruence of classical trajectories. The new equation we find is a 
Schrodinger-like one describing the evolution of the proper quantum state x- Neglecting 
higher order correction terms in fi, it reads 

-2ia'S'daX = H,x, H, = pi + pi + ^V{f3±). (4) 

Such an equation is in agreement with the assumption that the anisotropics describe a 
quantum subsystem of the whole Universe, i.e. that the wave function x depends only on 
f3± (in the Born-Oppenheimer sense). As matter of fact, the smallness of such a quantum 
subsystem can be formulated requiring that its Hamiltonian Hq is of order 0(e~^), where 
e is a small parameter proportional to fi. Since the action of the semi-classical Hamiltonian 
operator Ho = ci^da + a^U/n on the wave function •f' is of order 0{cr^), the idea that 
the anisotropics do not influence the isotropic expansion of the Universe can be formulated 
as Hq^/Ho^ = 0{e). Such a requirement is physically reasonable since, the semi-classical 
proprieties of the Universe as well as the smallness of the quantum subsystem, are both 
related to the fact that the Universe is large enough [2]. 

A purely Schrodinger equation for the wave function x is obtained taking into account 
the tangent vector to the classical path. Using pa ~ S' , the equations of motion (3) and 
considering the time gauge da/dt = 1, is possible to define the new time variable r such 
that dr = {NK/a^)da. In the asymptotic interesting region {a^ Ia = lj\/A) the evolution 
equation (4) rewrites as 

id^X = [-^0 + ^^(/3±)) X, (5) 

where r = (K/12\/^)a -|- 0{a^''). This is the Schrodinger equation for the wave function 
X describing the quantum variables /3±. The wave function (2) defines a probability distri- 
bution p{a,f3±) which appears to be p{a,f3±) = po{a)px{a, f3±), where po{a) is the classical 
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probability distribution for the semi-classical variable a and = |x|^ denotes the proba- 
bility distribution for the quantum variables /3± on the classical trajectories (3) where the 
wave function x can be normalized. 

In order to enforce the idea that the anisotropics can be considered as the only quantum 
degrees of freedom of the Universe, we consider the quasi- isotropic regime \P± \ 1. More- 
over, since we are interested at the link between the isotropic and anisotropic dynamics, the 
Universe has to be get through to such a quasi-isotropie era. In this regime, the potential 
term reads V{P±) = 8(/3+ + P-) + 0(/9^) and the Sehrodinger equation (5) can be then 
written as 

id^X = I + u>\r){pl + (3'_)) X, (6) 

where a;^(r) — C jr'^^'"' and C being a constant given by 2C — 1/6'*'^'' (k/I)^'''*. In other words, 
we are dealing with a time-dependent bi-dimensional harmonic oscillator with frequency 
a'(r). The quantum theory of an harmonic oscillator with time-dependent frequency is well 
known [9] and the solution of the Sehrodinger equation (6) can be analytically obtained. 
Through the introduction of the generalized invariant state, whose eigenstates are connected 
with those of a time-independent harmonic oscillator, and via an unitary transformation, 
the wave function Xn = X+X- reads 



X±=Xn{l3±,T)=A^^^hniP±/p)e^p J(pp"'+*p"')/3± • (7) 



In this formula A denotes the normalization constant, 7i„ are the usual Hermite polynomial 
of order n and p(r) and the phase Q;(r) are respectively given by 




)' «n = -(n+i 



P=l/^(1 + ^). ^n = -{n+^) I (8) 



It is immediate to verify that, when uj{t) ujo and p(r) ^ po = l/^/tJo (namely Q;(r) — > 
—ijJo{n + 1/2)t), the solution of a time-independent harmonic oscillator is recovered. 

Let us now investigate the probability density to find the quantum subsystem of the 
Universe at a given state. The anisotropics appear to be probabilistically suppressed as soon 
as the Universe expands enough far from the cosmological singularity (which we remember 
appears for o — » or t ^ oo). Such a feature can be immediately observed from the behavior 
of the squared modulus of the wave function (7) which is given by 

\Xnf ~ ^\hn+{p+/p)f\K_{p-/p)fe-^"/''\ (9) 

where = /3+ + /?- and with we omit the normalization constant. This probability 

density is still time-dependent through p = p(t) since the evolution of the wave function 
X is not traced by an unitary time operator. As we can see from (9), when a large enough 
isotropic cosmological region is considered (namely when the limit a cx) or r ^ is taken 
into account), the probability density to find the Universe is sharply peaked at the isotropic 
configuration, i.e. for \j3±\ ~ 0. In this limit (which corresponds to p ^ 0) the probability 

density |xn=o|^ of the ground state (n = r+ + n_ =0) is given by |xn=oP — 0), thus 
is proportional to the Dirac (5-distribution centered on /3 = (see Fig. 1). 

Summarizing, when the Universe moves away from the cosmological singularity, the prob- 
ability density to find it is asymptotically peaked around the closed FRW configuration. Near 
the initial singularity all values of anisotropy /3 are almost equally favored from a probabilis- 
tic point of view. On the other hand, as the radius of the Universe grows, the isotropic state 
become the most probable one. (For other similar approaches see [10].) 
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Fig. 1 The ground state of the wave function x{P± > ''') f^^r from the cosmological singularity, 
i.e. in the r — > limit. In the plot we take C = 1. 

It is worth noting that the key feature of such a result relies on the fact that the 
isotropic scalar factor a was considered as an intrensically different variable with respect 
to the anisotropics. It was treated as a semi-classical variable while only the two physical de- 
grees of freedom of the Universe (/3±) were described as real quantum coordinates. This way, 
a positive semidefinite probability density for the wave function of the quantum subsystem 
of the Universe can be constructed and a clear interpretation of it considered. The validity 
of such an assumption can be verified from the analysis of the Hamilton- Jacobi equations 
(3). In particular, the WKB function 'I'o — exp(iS' -I- In ^4) approaches the quasi-classical 
limit e*^ as soon as a ^ Ia {I a being the inflation characteristic length [7]). To corroborate 
the model, we have studied the classical limit too. Splitting the 5-function in two terms, 
respectively for the time variable a and for the anisotropies l3±, we achieve, if a 3> /yi, an 
analogous behavior of the anisotropies, i.e. them go to reduce themselves once one moves 
away from the cosmological singularity. 
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